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Neutron star atmosphere composition: the quiescent, 
low-mass X-ray binary in the globular cluster M28 
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ABSTRACT 

Using deep Chandra observations of the globular cluster M28, we study the quiescent 
X-ray emission of a neutron star in a low-mass X-ray binary in order to constrain the 
chemical composition of the neutron star atmosphere and the equation of state of dense 
matter. We fit the spectrum with different neutron star atmosphere models composed 
of hydrogen, helium or carbon. The parameter values obtained with the carbon model 
are unphysical and such a model can be ruled out. Hydrogen and helium models 
give realistic parameter values for a neutron star, and the derived mass and radius 
are clearly distinct depending on the composition of the atmosphere. The hydrogen 
model gives masses/radii consistent with the canonical values of 1.4 Mq and 10 km, 
and would allow for the presence of exotic matter inside neutron stars. On the other 
hand, the helium model provides solutions with higher masses/radii, consistent with 
the stiffest equations of state. Measurements of neutron star masses/radii by spectral 
fitting should consider the possibility of heavier element atmospheres, which produce 
larger masses/radii for the same data, unless the composition of the accretor is known 
independently. 

Key words: equation of state - stars: neutron - globular clusters: individual (M28 
or NGC 6626) - X-rays: binaries - X-rays: individual (CXOGlb J182432. 8-245208) 



1 INTRODUCTION 

Neutron stars (NS) are composed of the densest form of 
matter known to exist in our Universe, providing us with 
a unique laboratory to study cold matter at supra-nuclear 
density. For example, it is still not well understood whether 
exotic condensates occur in the NS core. The chemical com- 
position of the outer envelope is also uncertain, as well as 
the symmetry energy, the behavior of superfluidity among 
neutrons and protons, and the conductivity of the NS crust. 
Measuring the masses or radii of these objects can lead to 
useful constraints on the dense matter equation of state 
(EOS), and g ive insights of the composition of NSs (see 
lLattimerll201Q for a recent review). 

Masses have been determined very accurately for a few 
dozen NSs in binaries containi ng pulsars (for a recent compi- 
lation, see lLattimer fc Prakashll2010r ). Measured masses for 
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a variety of NS systems range from 1.0 to 2.5 Mq, with a 
canonical mass of 1.4 Mq. In particular, the high NS mass 
measured with high precision for the pulsar PSR J 1614-2230 
(1.97 ± 0.04 Mq) brings doubts on the presence of hyper - 
ons or free quarks in NS interiors (|Demorest et alJ 12010). 
though th ose solutions are not completely ruled out (e.g. 
IWeissenborn et al.|[2uTll : iMassot et al.ll2012l ). 

For systems in which the NS is hot enough to emit de- 
tectable X-rays from the surface, due to youth or accretion, 
the X-ray spectrum can b e used to constrain both the ra- 
dius and mass of the NS llLattimerl l2010l: lOzel et all l2010l : 
ISteiner et al.ll20ld : lGallowav fc Lampell2012l ). The mass and 
radius of isolated NS or ones in transient low-mass X-ray 
binaries (LMXBs) can be inferred from spectral modeling 
if their distances are accurately determined. In the case of 
accreting NS transients located in globular clusters (GCs), 
relativel y accurate distances are known (errors of ~5%, 
see e.g. iKrauss fc Chaboverj 120031 ). Even though, the un- 
certainties for individual NS measurements are still large 
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(iRutledge et all 12002 j; iGendre et al.l l2003al; JHeinke et all 



l2006l ; IWebb fc Barretll2007l ; iGuillot et al.ll201ll ). but the en- 
semble of observations can ult imately improve c onstraints 
on the dense matter EOS (e.g. ISteiner et alj|2010f ). 

It has been shown that the surface of a weakly mag- 
netic (B < 10 10 G) NS should be chemically very pure and 
dominated by the lightest element present as the heavier el- 
emen ts settle out of the atmosphere within second s to min- 
utes (JAlcock fc IllarionovllT980l : lBrown et al.ll2002l ). If there 
is accretion after the NS formation, the atmosphere could 
be composed of hydrogen -H- or helium -He- as heavier 
elements a re expected to be de s troyed via nuclear spalla tion 
reactions l|Bildsten et al.lll992l ; Ichang fc Bildstenl 120041 ). A 
fraction of the incident He also suffers spall ation reactions 
and may reform through fusion reactions IjBildsten et alJ 
Il993h . The ratio of H to He is thus not well determined. If 
no accretion takes place or if a ll lighter elements a re burned, 
heavy elements are expected IjChang et alJl20ld . and refer- 
ences therein). 

Different NS atmosphere models have been developed, 
but most recent work for low magnetic fields has focused on 
a pure H model, such a s the o nes developed by Zavlin et al.l 
l|l99rJ l. lGansicke et all (|2002l ), or lHeinke et"atl I2OO6). The 
latter model, NSATMOS, was further developed to repre- 
sent atmospheres of pur e He, carbon, nitrogen, oxygen or 
iron (|Ho fc Hei nkc 2009). In particular, such models were 
used for the low magnetic field NS located at the center of 
the Cassiopeia A supernova remnant, which wa s shown to 
harbor a carbon atmosphere (|Ho fc Hein ke 2009). 

The GC M28 (NGC 6 626) is located at a distance of 
D = 5.5 ± 0.3 kpc (from lHarris! Il99rj . l201p| . using mea- 
surements in lTesta et all 1200 ll ) at RA = 18 h 24 m 32.81 s and 
Dec = -24°52'11.2" (J20 00). The reddenin g toward M28 is 
E(B -V) = 0.42 ± 0.02 l|Testa et all 1200 ll ). implying a H 
column density of AT H = (2.33 ± 0.12) x 10 21 cm 2 (using 
IPredeh l & Schmitt 1995 for conversion). 

iBecker et alJ (J2003J) have previously reported on a set 
of ~40 ks Chandra X-ray Observatory ACIS-S observations 
of M28 (Obslds 2683, 2684, 2685). They suggested that the 
luminous, soft Chandra source numbered 26 in their work 
(IAU-approved source name CXOGlb J182432.8-245208) 
is a transiently accreting NS in a LMXB in quiescence 
(qLMXB). We keep the name source 26 throughout the 
text. F itting its spectrum with a H atmosphere model (NSA 
model, IZavlin et al.ll 19961 ) yielded the effective temperature 
Too ~ 90 eV and the projected radius Roo ~ 14.5 km with 
the mass set to 1.4 M@, corresponding to a NS temperature 
Tns ~ 125 eV and radius Ryis ~ 10.4 km. 

In this paper we focus on this NS qLMXB in M28 using 
one of the deepest Chandra observations of a GC to date and 
recent atmosphere models with different chemical composi- 
tions. We present the data in Section[2l and a variability and 
spectral analysis in Section [3] We finally discuss the results 
and implications in Section U 



2 DATA 

Th e Chandra data set u sed here was already presented 
bv lBogdanov et al.l 1120111 ). Two long observations were ac- 
quired on 2008 August 7 (Obsld 9132) and 2008 August 
10 (Obsld 9133) for 144 and 55 ks, respectively (PI: J. 



Gr indlay). We al s o inclu ded the three observations presented 
by IBecker et all (|2003l , Obslds 2683, 2684, and 2685, PI: 
W. Becker). All observations were taken with ACIS-S in 
VFAINT telemetry mode, and TIMED read mode with a 
frame time of 3.1 s and reading time of 0.041 s. The data 
re-processing, reduction, and analysis were perfor med in the 
same way as presented bv lBogdanov et al.1 i|201ll ). but using 
the more recent CIAO 4.4 and CALDB 4.4.8. In particular, 
we did not used the back-ground cleaning algorithm specific 
to the VFAINT telemetry mode as this procedure tends to 
reject real source counts for relatively bright sources. 

We us ed the CIAO tool wavdetect and the script 
PWDetect IjDamiani et al.lll997l ) to detect source candidates 
in the field of view. Subsequently, w e employed the IDL 
script ACIS Extract (|Broos et al.ll2010h to confirm the valid- 
ity of the source detections and refine the source positions. 
This lead to the detection of 101 sources inside the half mass 
radius of M28, and 35 inside the core radius, with a mini- 
mum of 4 counts per source. 

The qLMXB candidate, source 26 of IBecker et al.l 

(2003), is detected with a total of 10332 counts 



(~0.043 cts s _1 ) in the 0.3-6 keV energy band. We extracted 
the events of the target in a 1.25" radius region, enclosing 
95% of the total source energy for an effective energy of 
1 keV. We estimated that the contamination from close- 
by detected sources is <1%. The background was extracted 
from three 24" radius regions surrounding the GC core ra- 
dius, making sure the regions are source free. The data is 
clean from background flares, except in Obsld 2683 where 
the background flux reaches 5 times its median value for 
4 ks over the 14 ks observation. However, this corresponds 
to a negligible contamination of ~0.5 count in the aperture 
of the target (0.1% of the target counts in that Obsld), we 
thus kept all the data in this obsid. 

We used dmextract to extract a lightcurve with a 1000 s 
binning (minimum of 30 counts per bin), and dmcopy to ex- 
tract an event list for the source and background. Using 
specextract, which includes an energy dependent aperture 
correction (arfcorr), we extracted a spectrum and generated 
response files for each dataset. The spectra were binned with 
a minimum of 50 counts to use x 2 statistics for the fitting. 
We discarded bins flagged as bad, and bins below 0.3 keV 
where the response of the instrument is not calibrated (cal- 
ibration down to 0.3 keV with ACIS-S3 has been clearly 
improved in CALDB 4.43). 



3 RESULTS 

3.1 Light curve and search for variability 

We show in Figure [l] the lightcurve of the target for all five 
observations with a 1000 s binning, as well as the median 
energy of the source in each bin. We report the lightcurve 
characteristics and variability tests in Table [l] For each 
lightcurve we computed the \ 2 of a fit with a constant and 
performed a Kolmogorov-Smirnov (KS) test of variability. 
The source showed no significant variability in all Chan- 
dra observations. There is no feature in the light-curve that 
would indicate a periodic variability, or flares, and we found 
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Figure 1. X-ray lightcurve of the M28 Chandra source 26. The five Obslds are presented in consecutive pannels. The count rate and 
the median energy are calculated in the 0.3-6 keV energy band. 



Table 1. Properties of the X-ray observations of the M28 Chandra source 26. 



Obsld Date Exp. Time Counts 

[ks] 



Rate 

[icr 2 c/s] 



Median Energy 
[keV] 



Mean Energy 
[keV] 



X 2 /dof 



KS prob. 

% 



Frac. rms 

% 



2684 


2002-07-04 


12.2 


561 


4.4±0.2 


1.11±0.03 


1.2±0.5 


10.7/13 


22 


<23 


2685 


2002-08-04 


13.4 


569 


4.3±0.2 


1.08±0.02 


1.2±0.5 


12.8/14 


32 


<24 


2683 


2002-09-09 


14.3 


684 


4.9±0.2 


1.10±0.03 


1.2±0.5 


12.8/14 


48 


<21 


9132 


2008-08-07 


144.2 


6182 


4.29±0.05 


1.13±0.01 


1.2±0.5 


120.6/144 


52 


<12 


9133 


2008-08-10 


54.7 


2336 


4.3±0.1 


1.14±0.01 


1.2±0.5 


26.0/56 


15 


<21 



All 



238.8 



10332 



4.29±0.04 



1.13±0.01 



1.2±0.5 



173.3/200 



38 



<11 



no correlation between the flux and median energy of the 
source as a function of time. We also give the fractional rms 
(root mean square) of the light-curve. As the observed vari- 
ability is consistent with noise, a 3a limit on the fractional 
rms is given in Table [l] 

The possible 4<r increase of flux reported for Obsld 2683 
compared to 2684 and 2685 (|Becker et al.l 120031 ) is not ob- 
served in this work. The mean count rate is slightly higher 
(about 2a) in observation 2683, but the scatter around the 
mean for this Obsld is not unusual, as can be seen in Fig- 
ure [l] Similarly, we found no significant change in the mean 
energy. There may be an increase of the median energy be- 
tween 2002 and 2008 of ~0.03 keV, which is marginal (less 
than 2a). 



3.2 Spectral fitting 

We fitted simultaneously the five spectra extracted from the 
five different epochs with Xspec 12.7.0e (jAr naud 199 6]), us- 
ing t he pure H atmosphere model NSATMOS ( Heinke et al.l 
l2006f ) and a photoelectric absorption Ah along the line o f 
sight (TBABS, with abundances from I Wilms et all [2000). 
We first tied all parameters between the different datasets, 
which seems a reasonable assumption as the source did 
not show significant variability between the different epochs 
(section [3T}. We fixed the distance to 5.5 kpc and the nor- 
malization to 1 (i.e. we assume that all the NS surface is 
emitting). We also kept at first the mass and radius to the 
canonical values, 1.4 Mq and 10 km respectively. This re- 



turned a marginally acceptable fit with a reduced \ {Xv) 
of 1.27 and 144 degrees of freedom (dof). Data bins at en- 



ergies higher than 2 keV show systematic positive residuals, 
mainly visible for the Obsld 9132 spectrum (see Figure [2]). 

Given the count rate and softness of the source, which is 
almost on-axis, we expect some level of pile-up in the spec- 
trum (5% using PIMMS 4.fl). This effect occurs when two 
or more photon events overlap in a single detector frame and 
are being read as a single event, creating a high energy tail in 
a contaminated spectrum. We used the pile-up model com- 
ponent available in Xspec (|Davisll200ll ) with a frame time 
set to 3.1 s and a free a parameter (related to the proba- 
bility of events being retained as a good grade after filter- 
ing). This model gave the best fit {xl/doi = 0.88/143), with 
a = 0.41 ± 0.15. We refitted the spectrum with the mass and 
radius parameters free to vary. The best fit model (x^/dof 
= 0.87/141) is obtained for A H = (2.5 ± 0.3) x 10 21 cm 2 , a 
temperature kT cS = 125 ± 40 eV, a mass M = 1.4ig;| M 
and a radius R — 9 ± 3 km. Errors are at 90% signifi- 
cance and we considered only masses higher than 0.5 Mq 
and radii higher than 6 km. The 0.3-6 keV absorbed flux 
of the source (after removing the pile-up effect) is then 
(1.8 ±0.2) x 10~ 13 erg s" 1 cm -2 , and the unabsorbed lu- 
minosity ~ 1.6 x 10 33 erg s _1 (at 5.5 kpc). This spectrum is 
reported in Figure [2] We then ran the command steppar and 
obtained confidence contours for the mass and radius of the 
NS, which are more instructive than the best fit parameter 
values and errors (see Figure [3l left). 

Though no obvious variability was found in section l3Tl 
we tested this assumption by untying the normalization for 
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Figure 2. Spectrum of the M28 Chandra source 26 fitted with 
a NS atmosphere. The best fit model for the pure H atmosphere 
is reported as a dashed line and the dotted line shows the model 
without pile-up. Residuals are shown for the initial H atmosphere 
fit without pile-up and for different chemical compositions. 



all five Obslds, and only froze the Obsld 9132 value to 1. 
This led to the following best fit values for the normalization 
of each Obsld with their 90% confidence errors: 0.93 ± 0.08 
(2684), 0.89 ± 0.08 (2685), 0.99 ± 0.08 (2683), TO (9132), 
0.99 ±0.05 (9133). The largest deviation is thus at 2.3 signia 
(Obsld 2685), still consistent with unity. We performed a 
similar test with Nu and found no significant variability as 
well. Consequently, we confirm that no significant variability 
is observed between the different epochs, both in the flux and 
in the spectral shape. 

We also tested for the presence of a hard tail in the X- 
ray spectrum by adding a power law component. However, 
the normalization of this component quickly tends to zero, 
and the fraction of flux in the 0.3-8 keV energy band of this 
component cannot be higher than 5% or 7% of the total 
model flux (at the 95% confidence level) for a photon-index 
of 1.5 and 2.5, respectively. 

We performed the same fitting procedure with an at- 
mosphere model composed of pure He (usin g opacity ta- 
bles computed by the Opacity Projeco see iHo fc Heinkd 
120091 for details), and including the pile-up model. A sim- 
ilar good fit was obtained (xt/doi = 0.88/142, see Fig- 
ure 0) with iV H = (2.65 ±0.25) x 10 21 cm 2 , a temperature 
feTeff = nOtfo eV, a mass M = 2.0±?;jj M Q and a radius 
R = 14JI S km. The confidence contours obtained with the 



step-par command are reported in Figure [3] (right). We note 
that the regions delimited by the contours are not consistent 
at the 80% confidence level with the contours obtained with 
the H model. 

Finally, w e performed a simil ar fit with a carbon atmo- 
sphere model (|Ho fc HeinkdbOOfjf ). We obtain an acceptable 
fit (xl/dot = 0.88/142, see Fig ure H) but the paramete r 
values are excluded by causality (jRhoades fc Rufnnilll974T ): 
M > 2.6 M© for 7? = 10 ± 2 km. 



4 DISCUSSION 

The lack of variability of the source over six years agrees 
with the idea that the tar get is a LMXB in a quies- 
cent, low-accretion state (e.g. iRutledge et al.l|2002bD . Such 
qLMXBs can be X-ray transients with episodes of higher 
accretion that would heat up the NS star sur f ace t o 
the current observed temperatures (I Brown et al.l |l998|). 
Such a be havior has been observed for similar sou r ces in 
GCs (e.g. IRutledge et~a?1 l2002al: iHeinke et aI1l2003l . l20ld : 
iPoolev fc Hutll200rj ; iDegenaar fc Wiinandsll201ll V It is thus 
possible that this source becomes suddenly more luminous 
in the future. 

For both a H and a He model, we found good fits with 
an absorption consistent with the expected absorption from 
the GC reddening, suggesting that the source is located 
in the core of M28 with no or very low intrinsic absorp- 
tion. JTJiejn^sjuid radius are as expected for a typical NS 
(e.g. [Lattimcr 2010), and the temperature is in the expected 
range for qLMXBs. No hard tail is observed in the spectrum, 
which seems to be a commo n feature of these obje cts in GCs 
iJHeinke et al.ll2003J . but see lWiinands et alJl200ai . 

The only striking difference is that H and He atmo- 
sphere models give distinct contour regions of masses/radii 
at the 80% confidence level (Figure [3}. On the one hand, 
the H model gives a mass and radius consistent with the 
canonical value of 1.4 Mq and 10 km, and allows for the 
presence of exotic matter inside NSs (hyperons, quarks). On 
the other hand, the He model provides solutions with higher 
masses/radii, consistent with the stiffest EOS for NS interi- 
ors, most of them composed of neutrons and protons. 

The acc urate measurement o f 1.97±0.04 M for PSR 
J1614-2230 |Pemorest et al.ll2010T ) does not, by itself, indi- 
cate that a H atmosphere is ruled out since there are EOS 
that have a maximum mass above 2 Mq and lie within the H 
model confidence conto urs (AP4, but al s o others not shown 
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Lattimerl l2010l : IWeissenborn et al.l 
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l201ll : iMassot et al.ll2012T l. The He atmosphere model also 
allows for several stiff EOS that reach masses higher than 
2 Mq. We note that there are additional sources of error to 
consider (beyond the statistical uncertainties from spectral 
fitting) that does not allow for strict constraints on the NS 
parameters and the EOS, and may slightly enlarge the re- 
gions showed in Figure [3] Uncertainties on the distance were 
not included in our analysis, but this value is relatively well 
constrained as the source is in a GC. This has no impact 
when comparing H and He models as the same distance is 
used for both models. Although pile-up is a small effect and 
was carefully modeled in the fit, the additional parameters 
of the pile-up model may introduce additional systematic 
uncertainties. 
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Figure 3. Confidence levels for the mass and radius of the M2 8 NS Chandra source 26, u sing a H (left) or He (right) atmosphere model. 
A representative selection of EOS are reported (labelled as in lLattimer fc Pr akash 2001). The parameter s were not allowed to var y in 
the area "Not tested with the atmosphere model". We report in dark gray the area excluded by causality (Rhoades fc Ruffini 1974). 



The composition of the NS atmosphere depends on the 
accreting material, physical processes occurring during the 
accretion, and conditions on the NS surface. A non-evolved 
star will produce mostly H, which will quickly stratify to 
provide a pure H atmosphere. White dwarf donors (in so- 
called ultra-compact LMXBs) will provide mostly He, C/O, 
or O/Ne/Mg depending on the white dwarf. Ultra-compact 
LMXBs are observed to be muc h more common in GCs 
than in the rest of the Galaxy (jDeutsch et al.l 120001 ). Of 
16 bright LMXBs in 13 clusters, we have 11 orbital pe- 
riod measu rements, of which 5 indicate ultra-compact sys- 
tems (e.g. IZurek et all 120091 : lAltamirano et al.l l20ld ). In 
the rest of the Galaxy, only 9 ultra-compact systems are 
known among the ~80 bright LM XBs with period measure- 
ments JRitter fc Kolbll2003l . [2oTol ). Dynamical formation of 
ultra-compact LMXBs in GC core s explains this difference 
i|Verbuntlll987l ; llvanova et al"]|2005D . 

It is uncl ear whether spallation always produces H dur- 
ing a ccretion (|Bildsten et al.lll99l Il993l : IChang fc Bildstenl 
120041 ) . Theoretical work is needed to clarify the conditions 
for spallation. Obtaining a high-quality X-ray spectrum 
of an neutron star ultra-compact LMXB in quiescence at 
known distance would help clarify this question. 

Among the qLMXBs in GCs that were used to de- 
rive constraints on the mass and radius of their NS us- 
ing a H atmosphere models, some were reported to have 
a low mass or radius. For example, the source in the core 
of NGC 6553 (XMMU J180916255425) gives a NS radius 
Ans = 6.3lo| km for a m ass of 1.4 Mm (iGuillot et al.1 
l201ll ). The aLMX B in M13 llGendre et al.N20o"3al ) and the 



one in NGC 2808 (|ServiHat et al.ll2008bl lal) also showed rel- 
atively low masses/radi i with a radius lower t han 10 km or 
a ma ss around 1 M (|Webb fc BarretJ 120071 : ISteiner et al.l 
l2010h . however the values are not well constrained. In a 
same way, the qLMXB U24 in NGC 6397 was reported 



to show a relativ ely small radius J?n s — o.e>_ 0l . 

mass of 1.4 M© (|Guillot et al.ll201ll ). On the other hand, 

X7 in 47 Tuc was reported with higher values of mass/radius 

(R > 12 km or M > 2 M^. lHeinke et al.ll2006l). while the NS 

in to Cen appears as an int ermediate case IjRutledge et al.l 

l2002al : iGendre et ai1l2003bl ). as is the NS studied in this 

work. 

Following our study of the qLMXB in M28, it is pos- 
sible that some of those sources harbor a NS with a He 
atmosphere, rather than a H atmosphere. This would favor 
higher radii and masses for NS, and thus stiffer EOS, in 
agreement with the precise meas urement of relatively high 
masses for some NS (e.g. ~2 Mm. lDemorest et al.ll20101 ). We 
will thus try fitting other quiescent LMXBs with He (and 
carbon) atmospheres in future work. 

Identifying the composition of the atmosphere of known 
quiescent LMXBs is clearly of key importance, and we sug- 
gest three means of doing so. i) Spectroscopy, or (less time- 
consuming) narrow-filter photometry of optical counterparts 
can identify Ha emission from LMXBs in quiescence or out- 
burst and thus the presence of H; the LMXB in ui Cen 
l|Haggard et al.ll2004l ) and X4 and X5 in 47 Tuc (van den 
Berg et al., in prep) therefore possess H atmospheres, ii) Or- 
bital periods differentiate between ultra-compact and longer 
period systems; we note that long periods are known for X5 
and W37 in 47 Tuc (Heinkc et al. 2005), suggesting a main- 
sequence companion and accretion of H. iii) Finally, ther- 
monuclear bursts can distinguish between H-rich and H-poor 
environments, particularly at low (< 0.01 MEdd) accretion 
rates where H should bur n unstably (e.g. iFuiimoto et al.l 
ll98ll : lGallowav et alJlJOOsT ). 

This last point is of particular interest for the M28 
qLM XB, since a peculiar X-ray b urst was observed from this 
GC (iGotthelf fc Kurkarnilll997l ). This burst was unusually 
low-luminosity, suggesting burning on only one patch of the 
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star. The short timescale of this burst (r=7.5 s) requires He 
burning without the presence of H, and thus (given the qui- 
escent state) pure He accretion and a pure He atmosphere. 
Unfortunately we cannot be certain that this burst origi- 
nated from the known qLMXB, as other qLMXBs may be 
hidden among the fainter sources in this cluster. 
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